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ABSTRACT 

Massive stars, possibly red supergiants, which retain extended hydrogen envelopes until core col- 
lapse, produce Type II Plateau (IIP) supernovae. The ejecta from these explosions shock the cir- 
cumstellar matter originating from the mass loss of the progenitor during the final phases of its life. 
This interaction accelerates particles to relativistic energies which then lose energy via synchrotron 
radiation in the shock-amplified magnetic fields and inverse Compton scattering against optical pho- 
tons from the supernova. These processes produce different signatures in the radio and X-ray part 
of the electromagnetic spectrum. Observed together, they allow us to break the degeneracy between 
shock acceleration and magnetic field amplification. In this work we use X-rays observations from the 
Chandra and radio observations from the ATCA to study the relative importance of processes which 
accelerate particles and those which amplify magnetic fields in producing the non-thermal radiation 
from SN 2011ja. We use radio observations to constrain the explosion date. Multiple Chandra obser- 
vations allow us to probe the history of variable mass loss from the progenitor. The ejecta expands 
into a low density bubble followed by interaction with a higher density wind from a red supergiant 
consistent with Mzams ^ 16Mq. Our results suggest that a fraction of type IIP supernovae may 
interact with circumstellar media set up by non-steady winds. 

Subject headings: Stars: Mass Loss — Supernovae: Individual: SN 2011ja — shock waves — circum- 
stellar matter — radio continuum: general — X-rays: general 
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1. INTRODUCTION 

Type IIP supcrnovae display prominent P Cygni fea- 
tures arounf the time of peak luminosity, produced by 
hydrogen lines and their optical light curve plateaus 
for ~ 100 days in the rest-frame of th e supernova 
(jPoggett fc BranchlHOSllArcavi et al.ll2012[ ). This char- 
acteristic phase in their optical light curves is attributed 
to their progenitors retaining extended hydrogen en - 
velopcs until the time of core collapse. IPopovl (|1993f ) 
found that the duration of the plateau phase has a strong 
dependence on the mass of the hydrogen envelope and 
weak dependence on the explosion energy and the initial 
radius. These lines of evide nce and direct pre-explosion 
imaging (jSmartt et al.|[2009| ) suggest that these s tars ex- 
ploded as red supergiants. iSmartt et al.l ()2009[ ) found 
that two-thirds of the core collapse supernovae in their 
sample, volume lim ited to c? < 30 Mpc, are type IIP. 
iSmith et al.l ()201lD estimate the fraction to be closer to 
half. 

Red supergiants have been found insi de the Local 
Group with masses up to 25 Mq, but ISmartt et al.l 
(|2009( ) did not find any red supergiants with masses 
greater than 17Mq as progenitors of type IIP supernovae. 
Many solutions have been sugge s ted fo r the red super- 
giant problem. lO'Connor fc OttI (|2011f ) have suggested 
that the zero-age main sequence (ZAMS) mass, metallic- 
ity, rotation and mass-loss prescription controls the com- 
pactness of the stellar core at bounce which determines 
whether a core-collapse superno va will fail and instead 
form a stellar-mass black hole. IWalmswell fc Eldridgd 
(|2012D have suggested circumstellar dust as a solution to 
the problem of the missing massive progenitors. In this 
situation, understanding the nature, amount and vari- 
ability of mass loss from the progenitors of type IIP su- 
pernovae is crucial for resolving this puzzle. 

The ejecta from these explosions shocks the circum- 
stellar matter set up by the mass loss of the progenitor 
during the final phases of its life. Since the ejecta (~ 10^ 
km s~^) moves about a thousand times faster than the 
stellar wind (~ 10 km s~^), this expanding ejecta probes 
a millennium of red supergiant mass loss history in a 
year, a timescale which would otherwise be inaccessible 
in human lifetimes. This interaction accelerates particles 
to relativistic energies, which then lose energy via syn- 
chrotron radiation in the shock-amplified magnetic fields 
and inverse Compton scattering against op tical photons 
from the supernova. I Chevalier et al.l (|2006D have shown 
that these processes produce separate signatures in the 
radio and X-ray part of the electromagnetic spectrum. 
ICliakrab orti et al.l (|2012( ) have demonstrated that com- 
bining radio and X-ray spectra allows one to break the 
degeneracy between the efficiencies of shock acceleration 
and field amplification. In this work we use X-rays ob- 
servations from Chandra and radio observations from the 
Australia Telescope Compact Array (ATCA) to study 
the relative importance of particle acceleration and mag- 
netic field amplification f or producing the non therma l 
radiation from SN 201 Ija. iDwarkadas fc Gruszkol ()2012l ) 
have indicated that the expansion and density structure 
of the circumstellar matter must be investigated before 
assumptions can be made of steady wind expansion. It 
has been shown that the X-ray observations of SN 2004dj 
suggest variable mass loss though they do not rule out 



TABLE 1 

Observation of SN 2011,ia with Chandra 



Date 


XB Flux (0.3-10 keV) 
(10^^* orgs cm~^ 


SN Flux (0.3-10 kcV) 

(^IQ-li CTgg cm~^ S"-'-) 


2000 Jan 27 


1.01 ±0.11 


none 


2012 Jan 10 


0.81 ±0.10 


0.98 ± 0.17 


2012 Apr 03 


1.01 ±0.11 


4.08 ± 0.42 



Note. — Fluxes are model dependent. X-ray binary is modelled 
as tbabs(diskbb) and su pern ova is modeled as tbabs(powcrlaw) in 
XSPEC. See subsection 13.11 for details. Fluxes reported in this 
table are from the full model and not corrected for absorption. 



a constant mass-loss scenario (jChakraborti et al.|[20Tl) . 
In this work we use a multiple Chandra observation of 
SN 201 Ija to establish variable mass loss from the pro- 
genitor. 

2. OBSERVATIONS OF SN 2011JA 

SN 2011ja occurred in the ne arby galaxy NGC 494 5 
at a distance of 3.36 ± 0.09 Mpc (iMouhcine et al.[l2005D 
The supernova was first reported in lMonard et al.l ( 20111 ) 
where it was noted that Monard observed the SN at 
14.0 magnitude (unfiltered CCD) on December 18.1 UT 
and Milisavljevic obtained a spectrum on December 19.1 
UT that matched the type IIP SN 2004e t about a week 
after maximum light. For SN 2004et iCrockett et"all 
(|2011t ) have used the difference between the pre and post- 
explosion, ground-based o bservations to deduce a pro- 
genitor mass of ^ 8 Mq. iJerkstrand et al.l (|2012( ) have 
found a progenitor mass of ~ 15 for SN 2004et from 
late-time spectral modeling while iSahu et "aD ()2006l ) had 
found ^ 15 Mq from light-curve modelling. We com- 
menced our multi-wavelength campaign following the 
discovery of SN 2011ja. Our observations in X-rays and 
radio, reported and used in this work, are described in 
detail below. 

2.1. Chandra X-ray Observations 

SN 201 Ija was observed by a Target of Opportunity 
proposal (PI: Ray, Cycle: 13, ObsID: 13791) on 2012 Jan- 
uary 10 and subsequently using Director's Discretionary 
Time (ObsID: 14412) on 2012 April 03 from the Chan- 
dra X-ray Observatory, were used on both occasions, 
without any grating, for 40 ks each. The supernova was 
clearly detected in both of these observations. We also 
analyzed a pre-explosion 50 ks observation of the field 
(PI: Madejski, Cycle: 1, ObsID: 864) to look for possible 
contamination. Details of our observations are listed in 
Table [B 

Before spatial and spectral analysis, we processed the 
data from different epochs separately but identically. We 
followed the prescription from the Chandra Science Cen- 
ter using CIAO 4.4 with CALDB 4.4.8. We fihered the 
level 2 events in energy to only select ones between 0.3 
keV and 10 keV. The selected events were projected on 
the sky and the region of interest (a box with 20" sides, 
centered on the optical supernova position) was identi- 
fied. We masked the region of interest and generated a 
light curve from the remaining counts. We used this 
background light curve to identify fiaring and further 
masked time-ranges where the background count rate 
was greater than 3 times the rms. This left us with a 
table of good time intervals which was then used to se- 
lect the reliable events. The steps followed up to here are 
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Total counts pixel"' 

Fig. 1. — Photon counts obtained in the last Chandra epoch. The 
filtered Chandra events have been binned into image pixels which 
have 0.25 ACIS pixel sides. Note the spectral extraction regions 
for SN 2011ja (lower left in red) and contaminating X-ray binary 
(upper right in blue). 

the same as followed in [Chakraborti et all (|20T1) . The 
pre-explosion observation revealed the presence of a con- 
taminating source 1.35" from the supernova. 

The spectrum of the contaminating source was ex- 
tracted and it seems to be an X-ray binary. We shall 
explore the nature of this source in an upcoming work. 
However, the contaminant could not be precisely local- 
ized in the pre-SN exposure as it was 4.6' away from 
Chandra's bore-sight, where the PSF degradation is sub- 
stantial. Thus, detailed analysis of the post explosion 
observations were required in order to separate the su- 
pernova's flux from that of the contamination. We fit, 
using Sherpa, the 2-D image created from the event file 
of each observation by binning over the region of interest 
in 0.25 ACIS pixel sizes (See Figure[T]). We also created a 
PSF image file to use as a template using the CIAO tool 
mkpsf, which extracts a PSF model image from a given 
standard PSF library hypercube given an energy, offset 
and sky or detector physical coordinates. A source model 
of two point sources and a fixed background, convolved 
with the appropriate PSF was then fitted to the image. 
This allowed us to determine the relative positions of 
the supernova and the contaminating source (just 1.35" 
away) which are found with sub-pixel accuracy in the 
image plane. The extraction regions were defined as cir- 
cles with radius 0.67" centered at these positions. Given 
that these sources are on axis in the ToO observation, 
this area should contain ^ 90% of the fiux. However, we 
have applied an energy-dependent point-source aperture 
correction, using the CIAO tool arfcorr, to account for 
any missing flux. Wc generated the spectra, response and 
background files separately for both sources. We did not 
bin the data over energy and instead used unbinned data 
for further analysis to fully exploit the spectral resolution 
of the instrument. 



TABLE 2 
ATCA Observations SN 2011ja 



Date 


Frequency 


Flux Density 




(GHz) 


(mjy) 


2011 Dec 19 


18.0 


0.53 ±0.09 


2011 Dec 20 


9.0 


0.85 ±0.11 


2011 Dec 20 


5.5 


0.54 ±0.10 


2012 Apr 11 


9.0 


< 2.1 



2.2. ATCA Radio Observations 

SN 201 Ija was observed in t he radio so o n af- 
ter discovery with the ATCA (jRvder et all 120 lit ) 
and the Giant Metrewave Radio Telescope (GMRT) 
(jYadav fc Chakrabortil |20T2| ). Table H lists the flux 
densities observed (or upper limit) with the ATCA us- 
ing the Cornpact Array Broad-band Backend (CABB; 
I Wilson et all (20 111 ) which provides 2x2 GHz IF bands. 
Total time on-source ranged from 1 to 2 hours, yielding 
sufficient uv-coverage to comfortably separate SN 201 Ija 
from the side-lobes of the radio-bright nucleus of NGC 
4945 some 250" to the northeast. The ATCA primary 
flux calibrator, PKS B1934-638 has been observed once 
per run at each frequency to set the flux scale at all fre- 
quencies. It also deflned the bandpass calibration in each 
band, except for 18 GHz where the brighter source PKS 
B1253-055 was used instead. Frequent observations of 
the nearby source PKS B1320-446 allowed us to monitor 
and correct for variations in gain and phase during each 
run, and to update the antenna pointing model at 18 
GHz. The data were edited and cali brated u sing stan- 
dard tasks in the MIRIAD package (jSault et a l. 1995^, 
and images made using robust weighting. Fluxes in Ta- 
ble [2] were derived using the uvfit task to minimize uncer- 
tainties introduced by cleaning, phase stability, etc while 
fltting in the image plane, an d the uncerta i nties calcu- 
lated in the same manner as IWeiler et al.l (|20Tll) . An 
upper limit of 3 mJy (3 cr) was obtained using GMRT 
observation on 2012 January 11 UT at an effective fre- 
quency of 1264 MHz. 

3. NON-THERMAL EMISSION 

The fast moving supernova ejecta shocks the slowly 
moving pre-explosion circumstella r matter set up by 
the stellar wind of the progenitor (iCheval erlll982[ ). In 
a type IIP supernova, (jChakraborti et aL " l2012| ) point 
out that the post forward shock circumstellar matter 
is at too high temperature and low density to pro- 
duce a significant thermal contribution to the Chan- 
dra fiux. However non-thermal electrons accelerated 
at the forward shock can produce mo st of the radio 
emiss ion seen in type II P supernovae (iChevalier et ahl 
I2OO60 . Either thermal dSr itaria ct al.| 120031 ) or non- 
thermal (jBiornsson fc Franss on 2004i)electrons can In- 
verse Compton scatter a fraction of the optical supernova 
photons into the Chandra X-ray band. A non-thermal 
electron population specified by an index p produces syn- 
chrotron emission in radio with spectral index [p — l)/2 
and inverse Compton scattered X-rays with photon index 
{p+l)/2. 

3.1. X-ray Spectral Fitting 

The supernova spectra were imported into XSPEC 
12.7.1 for spectral analysis. Spectra from both epochs 
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4x10^ 6x10^ 8x10^ 

Thermal Model Emission (APEC Normalization) 

Fig. 2. — Probability density contours for models fitted to the 
Chandra data. Thermal flux (APEC normalization) on x-axis and 
inverse Compton flux (power law normalization) on te y-axis. The 
thermal and inverse Compton fluxes are anti-correlated as their 
sum has to aecount for the combined flux observed with Chandra. 
Note however that the data does not rule out a zero thermal flux. 

were fitted simultaneously with the sum of a non-thermal 
inverse Compton component using powerlaw and ther- 
mal emission fr om a coUisionally-io nized diffuse gas using 
APEC models (jSmith et al.ll2001[ ). The spectra was at- 
tenuated by a Tuebi ngen-Boulder ( tbab s i n XSPEC) ISM 
absorption model (jWilms et all 120001) . IMorival (|2012f ) 
has shown that if the circumstellar medium is dense 
enough, coUisional ionization equilibrium can be estab- 
lished in the early stage of the evolution of the super- 
nova remnant, especially in the reverse shocked plasma 
(jChakraborti et al.|[20Ta) . The APEC plasma temper- 
ature was fixed at 1 keV as suggested by the tem- 
perature of th e reve rse-shocked material calculated by 
iNvmark et al.l (|2006l ) and demonstr ated to be appro- 
riate for type IIP supernovae by iChakraborti et all 
2012f ). The powerlaw photon index representing the in- 
verse Co mpton componen t was fixed at 2, as observed by 
ICliakrabo rti et al.l ()2012| ) corresponding to p = 3. The 
column density for absorption was kept free but pegged 
to be same at both epochs. The APEC emission mea- 
sure and powerlaw normalization were solved for both 
epochs from this joint analysis. Each spectral channel 
would have too few photons for an useful estimation 
because the fits were performed on unbinned data. So 
we used the (jUash I97?i) statistic to perform our fits. 

In order to determine if the free parameters were in- 
deed well constrained by the data, we ran Markov chain 
Monte Carlo (MCMC) simulations with 10000 steps 
over the multidimensional space of all the free param- 
eters. We show in Fig [2] that while the inverse Comp- 
ton component is well-detected, there is no conclusive 
evidence for the therm al plasma component . It has 
already been predicted (IChevalier et alj [20061 ) and ob- 
served (jChakraborti et al.ll2012[ ) that in tvpe IIP super- 
novae the early X-ray emission is dominated by the in- 
verse Compton component. Therefore for simplicity, we 
set the thermal flux to zero in the subsequent analysis. 
We generated 10000 simulated spectra for the best fit 
model to test its goodness of fit. ~ 50% of these spectra 
were found to have cstat less than that for the real data. 
This leads us to conclude that our best fit model provides 
a good fit to the data. The absorption column density 




Fig. 3. — X-ray spectra of SN 2011ja. Bars are counts from 
Chandra, binned for display. Dotted line is the powerlaw model 
for the inverse Compton flux. Solid line is the full model after 
absorption, and binning. Black is first epoch (January 10) and red 
is second epoch (April 03). Note the significant increase in flux for 
the second epoch. 

for the best fit model is iVn = (7.5 ± 0.9) x lO^i cm^^. 
Refer to Table [1] for the X-Ray fiuxes determined from 
these models at each epoch and Fig |3] for the spectra. 

3.2. Radio Spectral Fitting 

Radio emission from supernovae can be modeled 
as synchrotron emission from the interactio n between 
super nova ejecta and circumstellar matter (jChevalieil 
Il982| ) . Most radio light-curves show a power- law decline 
at late ti mes and rise due t o low-frequency absorption 
processes (jWeiler et al.ll2003) at early times. The rising 
part has been often modeled as free-free absorption al- 
lowing one to estimate the circumstellar density. On the 
other hand assuming synchrotron self-absorption yields 
an approximate radius of the emiss ion region at the time 
of peak fiu x both for Newtonian ([Chevalierl 119981) and 
relativistic (jChakraborti &: Ravll2011|) explosions. If an- 
other mechanism such as free-fee absorption is dominant , 
the radius must be even larger. IChevalier et~all ()2006l ) 
argue that the expansion velocities of ~ 10^ km im- 
plied for the Type IIP supernovae in the synchrotron self- 
absorption model are similar to those expected for cir- 
cumstellar interaction and are thus consistent with this 
absorption mechanism. 

We therefore fit the radio spectrum of SN 201 Ija with 
a synchrotron self-absorption model (see Fig H}. The 
spectral indices of the optically thick and thin parts are 
fixed to —1 and 5/2 respectively. The best fit determines 
the two free parameters, namely the peak fiux density 
{Fp = 0.829 ± 0.033 mJy) and the peak frequency {ly = 
9.29 ±0.39 GHz) of the spectrum on 2011 December 19. 

4. BLASTWAVE PARAMETERS 

We can now use the results of the multi- wavelength ob- 
servation and analysis described above, to constrain the 
parameters of the supernova blastwave and its interac- 
tion with the circumstellar matter. 

4.1. Explosion date 

It is important to determine the explosion date of SN 
201 Ija, so that meaningful comparison with models for 
circumstellar interaction arc possible. Constraining the 
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Fig. 4. — Synchrotron self-absorption model fit to the SN 2011ja 
flux densities observed with the ATCA. Note the optically thin 
part to the right (~ i'"^) and the optically thick part to the left 

explosion dates of nearby core collapse supernovae is also 
important from the perspective of multi-messenger as- 
tronomy. For example neutrino detectors have limited 
direction sensitivity and one would need constraints on 
the explosion date to discuss the possibility of associating 
a few (say ~ 2) neutrinos with a nearby supernova. 

The forward shock accelerates electrons to relativistic 
energies. Synchrotron losses from these electrons p rodice 
the radio emission from supernovae. IChevalied ([l98l 
modelled the radio emission by assuming that a frac- 
tion Ee or of the thermal energy is used to accelerate 
electrons and amplify magnetic fields respectively. Us- 
ing this assump tion and a self similar blastwave solution 
iChevalieil ()1998f ) derived the radius of the radio emitting 
region as 



-1/19 / p N 9/19 



VMpcy VSGHz 



(1) 



where Fp is the peak flux at peak frequency v, the 
equipartition factor is defined as a = Ce/tB the ratio 
of relativistic electron energy density to magnetic en- 
ergy density and / is the fraction of the spherical vol- 
ume occupied by the radio-emitting region. It can be 
seen from the equation above that the estimated radius 
is insensitive to the assumption of equipartition. We 
therefore use our radio spectrum to estimate a size of 

Rs = (1.848 ± 0.085) x 10 ^^ cm. 

We further note that iMilisavlievic eTall (poTTh re- 
ported that the absorption minimum of the Ha line is 
found to be blue-shifted by about 11000 km s~^, which 
is usual for type IIP supernovae. Assuming a 10% uncer- 
tainty in the expansion velocity, we get an age of 19.4 ±2 
days at the time of the radio observations. Our best es- 
timate for the explosion date is therefore 2011 November 
30 UT (±2 days). This is consistent with the fact that 
the optical spectra taken on December 19.1 UT matches 
best with the spectrum of type IIP eve nt SN 2004et taken 
abou t a week after maximum light (jMilisavlievic et al.l 

[2oii|). 



4.2. Circumstellar interaction 

The forward shock accelerates electrons to relativistic 
velocities and amplifies magnetic fields, which are respon- 
sible for the radio emission from supernovae. The super- 
nova ejecta collides inelastically with the circumstellar 
matter. The external density is described by a power 
law profile p oc r"'', where s = 2 for a steady wind. We 
therefore have 

p - ^ - (2) 

where M and are th e mass los s rate and velocity of 
the wind respectively. IChevalieil (jl982| ) sets the nor- 
malization o f the circumstell ar density profile ss A = 
M/(47r?;w). IChevalieil (11981) assumes that a fraction 
Ce of the thermal energy is used to accelerate electrons 
while a fraction eg is used to amplify magnetic fields. 
Hence these microphysical parameters determine the ra- 
dio brightness of a supernova, whi ch is not a direct mea- 
sure o f the circumstellar density. I Chevalier fc FranssonI 
(|2006l) calculate that that radio emission can constrain 
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(3) 



at a time 10 x iio days after explosion. Here 
eB_i = es/O.l and = ^/(5 x 10" gem"!) 

is a non-dimensiona l ized proxy for A defined by 
iChevalier fc FranssonI (j2006') . Our radio spectrum de- 
termines = 1.979 ± 0.398. 

The electron population that emits radio synchrotron, 
also inverse Compton scatters optical photons into 
X-rays. This process dominantes the non-thermal 
part of the X-ray spectrum of type IIP sup e rnova e 
during the plateau phase (jChakraborti et akl 120121) . 
IChevalier fc FranssonI (|2006[ ) have shown that the inverse 
Compton flux at 1 keV produced by accelerated electrons 
with p = 3, is given by 



E 



dE 



8.8 X 10367„„„ 

/ Uox{t) 

\ 1042 ergs s- 



no 



(4) 



Here the smallest Lorentz factor for an accelerated elec- 
tron is 7min and V4 = 1.1 is the expansion velocity in 
units of 10"* km s^^ at 10 x tio days. 

During the first epoch (tio ~ 1.94) of Chandra ob- 
servations, we find the inverse Compton flux density to 
be (7.27 ± 1.50) x 10^^ ergs s"!. This gives us the left- 
hand-side of Equation [H We use the observed value of 
Si, = 1.979 as found using our radio observations and V4 
(~ 1.1) seen in optical spectra. This gives us 



In 



,a"/i^ X 



ibol(i) 



10"^^ ej-gs g-i 



1.55. 



(5) 



If the spectrum of accelerated electrons starts from those 
at rest we would have 7min = 1- Following the work by 
IChevalier fc FranssonI ()2006D in the case of SN 2002ap, 
we consider relativistic electrons with 7niin ~ 2.5 and 
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a bolometric luminosity of 10'^^ crgss"^ as is usual for 
the plateau phase of type IIP supernovae. This gives 
us a ~ 0.44. This is a direct test of the equipartition 
assumption which demonstrates that the electrons and 
magnetic fields are not far from equilibrium. 

We can now use the above result, Equation |3] defining 
Si, and a characteristic value of es = 0.1 to constrain the 
pre-explosion mass loss rate of the progenitor to be 



Mes-i = (2.7 ±0.5) x 10" 



10 km s 



Mq yr- 



(6) 



during the last century before explosion. This is smaller 
than what is expected for typical red supergiant pro- 
genitors (see Figured]) but similar to the mass loss rate 
seen for the progenitor of SN 2004dj (|Chakraborti et aH 
[20T1 . 



5. TEMPORAL VARIATION 



The X-ray spe ctra of ty p e IIP supernovae, such 
as SN 2004et (IRho et al.l [2001 and SN 2004dj 



- (IR _ - . 

(jChakraborti et alT 20121), so f ten an d fall in luminosity 
over time. iChakraborti et al.l ()2012D have shown that at 
initially, the Compton flux dominated the spectrum and 
produces a harder spectrum. But the numcr density of 
seed photons available for scattering decreases when the 
optical luminosity from the supernova falls rapidly at the 
conclusion of the plateau phase. Therefore at late-times 
the thermal emission from the reverse-shocked plasma 
dominates the spectrum and makes it softer. We ob- 
tained a second epoch of Chandra observations of SN 
2011ja to study its temporal variation. 

5.1. X-ray rise 

The inverse Compton X-ray flux varies in time due 
to the expansion of the blastwave and the change in 
the number density of seed photons. If the blastwave 
encounters circumstellar matter set up by the uniforrn . 
wind of the progenitor, IChevalier fc FranssonI (|2006f ) 
have demonstrated that the inverse Compton flux varies 

Therefore during the plateau phase of a type IIP super- 
nova when there is a nearly constant Lboi the inverse 
Compton flux should be oc t^^. Since the supernova is 
likely ^ 125 days old during the second epoch of Chandra 
observations, one would expect an X-ray flux reduced by 
at least a factor of ^ 6.4. Instead we find an increase in 
flux by a factor of ^ 4.2. This is inconsistent with the 
predictions assuming a circumstellar density oc r~^ set 
up by a steady wind. Unless micro-physical parameters 
such as the efficiency of electron acceleration Ee changed 
between the two epochs, this implies a variable mass loss 
rate for the progen itor. A similar rise was reported by 
iPoolev et alT(|2002[ ) for SN 1999em at around ~ 100 days 
after explosion, where the total flux nearly doubled from 
the previous observation, despite the continued decline of 
the high-energy X-rays. Therefore the spectra softened 
remarkably during the sudden rise in flux. For SN 2011ja, 
the Chandra X-ray source hardness ratio, calculated as 
{H — S)/{H + S) where S and H are the counts in the 
0.5 to 2 keV band and the 2 to 8 keV band respectively, 
changes from 0.08 ± 0.20 to 0.10 ± 0.10 which is consis- 
tent with no change between ~ 41 and ~ 84 days after 



Me 



B-l 



explosion. If the increase in X-ray flux resulted from in- 
creased circumstellar interaction, then the inverse Comp- 
ton component should scale as oc M while the thermal 
component should scale as oc M^, eventually overtaking 
the former for a mass loss rate of ^ 1O~^M0 yr~^. The 
hardness ratio should have changed if for example ther- 
mal emission from the reverse shocked plasma started 
to dominate or if the absorption column got reduced. 
Therefore the emission continues to be dominated by the 
non-thermal inverse Compton component. This is con- 
flrmcd by our exploration of the model parameter space 
following the MCMC method outlined in Section [3Tl 

5.2. Density enhancement 

The supernova ejecta moves a thousand times faster 
than typical red supergiant winds and probes three cen- 
turies of mass loss history during the period spanned by 
our observations. This potentiallly offers us a glimpse 
into the last stages of stellar evolution before core col- 
lapse. Inverse Compton flux scales directly with the num- 
ber of seed photon (roughly constant during the plateau), 
mass loss rate and inversely with time. Therefore we 
argue that one needs an enhanced density by a factor 
of ^ 27 to account for the increased flux in the second 
epoch. We find a mass loss rate of 

(^•^^^•^)^^«"(io^) 

However, this calculation is based on the scaling rela- 
tion given in Eq. JH which was d erived from the self 
similar solution of iChevalieii ()1982|) . assuming a steady 
wind. In Eq. 2] the emission is also a funtion of the shock 
velocity, so if the supernova slowed down significantly be- 
tween the two epochs, then it will increase the density 
enhancement required. Therefore, while the argument is 
physically well-grounded, the result can be incorect by a 
factor of a few. This should motivate investigations into 
the propagation of supernova blastwavcs into environ- 
ments with density jumps. One such situation could be 
the evolution of sup ernovae in c i rcums tellar wind-blown 
bubbles explored bv iDwarkadad ()2005[ ). 

6. DISCUSSION 

The results presented here tell us about the circum- 
stellar environments of the progenitors of type IIP super- 
novae and the equipartition assumption often invoked in 
supernovae shocks. These are explored in brief below. 

6.1. Equipartition in radio supernovae 

Radio emission from non-thermal sources can be ex- 
plained with various amounts of acc elerated electrons 
and amplified magnetic fields (jPachol czvk 197'c|). It is 
however difficult to const rain the relative contribution 
of each. iBurbidg^ (jl956[ ) demonstrated that assuming 
minimum energy in the emission region, implies that 
these energy densities are ap proximately in equiparti- 
tion. iScott fc ReadheadI (|1977[ ) observed that for sources 
having low-frequency spectral turnovers, the total energy 
is within a fac t or of a few of the equipartition energy. 
iKulkarni et"all (|1998| ) used the equipartition argument 
to estimate t he radius e x pansi on of the GRB associated 
SN 1998bw. iChevaheri (119981 ) noted that the inferred 
radii of synchrotron self absorbed sources are insensitive 
to the assumption of equipartion. Therefore independent 
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Fig. 5. — Zero Age Main Sequence Mass (ZAMS) and wind 
mass loss rate (during the last 100 years of stellar evolution) 
for MESA runs (magenta ) , an d theoretical line (blue) from 
INieuwenhuiizen fc de Jagerl 1119901 ) (with R = lO^iJ©) plotted for 
comparison. Shaded boxes are la confidence intervals for the 
mass loss rate observed in SN 2011ja (corresponding to eg = 
0.1) first epoch (r e d), se cond epoch (green) and SN 2004dj from 
IChakraborti et al.l II2012I ) (grey). There is disagreement between 
the predicted and observed mass loss rate of the SN 2011ja pro- 
genitor during the first epoch and the subsequent agreement during 
the second epoch. 

size measurements do not tightly constra in the equiparti- 
tion parameter a. IChandra et al.l (j2004D have suggested 
looking for a spectral break so that the magnetic field 
and the size of the radio-emitting region are determined 
through unrelated methods. 

Using SN 2004dj as a prototype IChakraborti et aTl 
(|2012i) have demonstrated that if the inverse Compton 
scattered supernova photons can be detected in the X- 
rays along with radio synchrotron emission, the param- 
eters of the system can be derived without t he he lp of 
the equipartion argument. iSoderberg et al.l ()2012D has 
explained the radio and X-ray properties of SN 2011dh 
using the synchrotron and inverse Co mpton mecha - 
nisms respectively with a a ^ 30 wh ile IMaedal (poTl) 
finds a < 1 and iHoresh et al.l (fMl find a ~ 1000. 
IBarniol Duran et al.l ( 2013[) have similarly shown that for 
sources where the Synchrotron-self-Compton component 
can be identified, the equipartion argument is not neces- 
sary. In this work, by comparing the radio synchrotron 
and X-ray inverse Compton flux densities, we have deter- 
mined the equipartition factor for SN 2011ja. Our result 
shows that the plasma in SN 2011ja is indeed close to 
equipartition. 

6.2. Progenitors of type IIP Supernovae 

Massive stars (M > 8Mq) evolve from a main sequence 
blue giant to a red supergiant and then explode as su- 
pernovae. Such a sequence of events is consistent with 
observations of most type IIP supernovae, their progen- 
itors and circumstellar interaction. Here we compare 
the observed mass loss rates of the progenitor of SN 
2011ia, with those of its expected progenitors. MESA 
(jPaxton et al.|[20Tl was used to evolve stars with masses 
between 11 to 19 M (7i, for a metallicity of z = Q.5Zq. 
(jPaxton et al.|[20TTI ) Sec 6.6 describe the mass loss pre- 
scription used in our simulations as the Dutch Scheme. 
This prescription turns on a RGB wind at the cor- 
rect burning stage. Changes in surface temperatures 
in different evolutionary stages are also taken into ac- 



durm g tne nXjO pna.se w nicn follows tne prescriptior 
from Ide Jager et al.l ()1988D . We evolved the stars un- 
til they reached a central density of 10^^ g cm~'^. The 
averaged the mass loss over the final hundred years to 
obtain our fiducial values. The supernova ejecta should 
encounter this circumstellar matter over the first months 
after explosion. 

The results of the above simulation are compared with 
the observations in Figure [SJ For this comparison, recall 
that with time circumstellar interaction probes deeper 
into the progenitor's mass loss history. Note that the 
circumstellar density observed in the first epoch, cor- 
responding to the final stage of mass loss before core 
collapse, is fa r below expec t ed val ues for red supergiant 
progenitors. IChugai et al.l (j2007[) also noted a similar 
disagreement for mass loss rates obtained from optical 
spectroscopy of S N 1999em and SN 2004dj . Similar val- 
ues are derived by IChakraborti et al.l ()2012l) for the mass 
loss rate of SN 2004dj from X-ray spectroscopy. Note 
however that by the second epoch of observations, corre- 
sponding to an earlier stage of the progenitor star's life, 
the mass loss is consistent with that from a red super- 
giant of M > 16Mq. Similarly, iKochanek et al.l (|2012D 
have used X-ray observations of SN 2012aw to infer a 
varying mass loss rate. 

The mass loss rate implied by Eq. [6] is very low for 
a slow moving red supergiant wind, but not for a fast 
moving wind of a blue giant. At the same mass loss rate, 
a faster wind sets up a lower circumstellar density. For 
example, the progenitor of SN 19 87A was identified to be 
a blue giant (jArnett et al.l 1 19891 and references therein). 
The circumstellar environment of a red supergiant is set 
up by its slow and dense stellar wind. However, if the 
star becomes a blue giant again before it undergoes core 
collapse, it will blow a hot low density bubble within 
the red supergiant wind. In such a situation, the su- 
pernova remnant will encounter little circumstellar mat- 
ter at early times but have stronger interaction subse- 
quently. This is indeed consistent with our observation 
and anal ysis of SN 201 Ija. Such an interaction was pre- 
dicted bvlChevalier fc LiaMl ([19891 ) for SN 1987A, while 
iLuo &: McCravl ()199H ) predkted a sharp rise in radio and 
X- ray luminosities. This rise w as subsequently observed 
b ylStavelev-Smith et al.l (|1992D . 

iBauer et al.l (|2008| ) have suggested that a similar situ- 
ation can e xplain the observations o f the exotic type Iln 
SN 1996cr. iDwarkadas et al.l (|2010[ ) studied the increas- 
ing X-ray flux from SN 1996cr, using a hydrodynami- 
cal model, computing non-equilibrium ionization spec- 
tra and light curves, then fitting them to observations, 
to fully understand the circumstellar env ironment. Cir- 
cumstellar interaction in this and other (IChandra et al 

mtype Iln have received recent attention. lOfek et al 
have observed of a remarkable mass-loss event de- 
tected 40 days prior to the explosion of the Type Iln su- 
pernova SN 2010mc. Radio and X-ray observations of SN 
2003bg (Soderberg ct al.' 20061), the ordinary type Ic SN 
2007gr (S oderberg et al., 20 , 1 0| ) and the broad-lined type 
Ic SN 2007bg (fSalas et al.ll20l'3l ) ha ve revealed the i r com - 
plex circumstellar environments. iWellons et al.l ()2012D 
have explored the unusual circumstellar environments for 
type Ibc Supernovae 2004cc, 2004dk and 2004gq. The 
peculiar SN 1987A is also interacting with a complex cir- 
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cumstellar environment (|Dewev et al.l 120121 ) . Compared 
to type Iln and type Ibc supernovae, few type IIP su- 
pernovae have had their environments stud ied in detail 
like Supernov ae 1999em (jPoolev et al.ll2002| ) and 2004ct 
isra ei~all[20nl . Here we suggest that the regular 
type IIP SN 2011ja is undergoing interaction with a com- 
plex circumstellar environment set up by a non-steady 
wind. 

7. CONCLUSIONS 

Radio and X-ray observations of SN 2011ja allow us to 
measure micro-physical parameters such as the ratio of 
energies which goes into accelerating electrons and am- 
plifying magnetic fields. It is deduced that in this case, 
the plasma is not far from equipartition which is often 
assumed in the study of supernova circumstellar inter- 
action. Radio observations have allowed us to constrain 
the date of explosion. Multiple epochs of Chandra ob- 
servations have allowed us to demonstrate that the su- 
pernova initially encountered a low density region, incon- 
sistent with the expected mass loss rate of a red super- 
giant progenitor. The fast moving cjccta subsequently 
catches up to the slowly moving wind from a possibly 
-^ZAMS ^ I6M0 red supcrgiant. This interaction with 
a lower density region followed by stronger circumstellar 
interaction is consistent with a SN 1987A-like blue giant 
progenitor with a fast wind for SN 2011ja and indicates 
that a fraction of type IIP supernovae may happen inside 
circumstellar bubbles blown by hot progenitors or with 
complex circumstellar environments set up by variable 
winds. If a ris e at ^ 100 days in X-ray fluxes as seen 
in SN 1999em (iPoolev et al.|[200l and SN 2011ja (this 
work) is more common than expected, multiple Chan- 
dra observations of young nearby type IIP supernovae 



are needed in the first year after explosion. In such a 
situation, for studies of c ircumstellar interac tion of type 
IIP supernova e, optical (iChugai et al.l [20071 ). radio and 
X-ray ( Cheval ier et al.|[2006D spectra at multiple epochs, 
are necessary to constrain the environments and progen- 
itors of type IIP supernovae. 
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